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1 Introduction 



Recent years have seen a substantial growth of interest in the calculation of excited baryon spectra 
in lattice QCD [1-11] that has largely been motivated by the present experimental programs at 
Jefferson Lab [12] and ELSA [13] (for the latest lattice reviews, see e.g. [14-17]). In general, 
the procedure of extracting resonances from lattice QCD data differs from the one used in the 
stable particle case, since resonances do not correspond to isolated energy levels of the total 
Hamiltonian. The standard approach, originally proposed by Luscher [18-21] (see also [22-24]), 
is based on studying the volume dependence of the spectrum being determined by the two-body 
scattering phase shift in the infinite volume. Near the resonance energy, where the phase shift 
rapidly passes through tt/2, an abrupt rearrangement of the energy levels known as "avoided level 
crossing" takes place. It has been argued that the observation of this phenomenon in lattice data 
can serve as a signal of the presence of a resonance and enables one to determine its parameters. 
In Refs. [25-27] the approach has been further generalized for moving frames. Note also that 
Liischer's approach has been recently applied to study nucleon-nucleon phase shifts at low-energy, 
as well as two-body shallow bound states [28-32]. 

As alternative approaches to this procedure we mention, e.g., Ref. [33], where it has been 
shown that the presence of a narrow excited state above the threshold modifies the simple ex- 
ponential decay law of the time-sliced two-point function. The decay width within this ap- 
proach is extracted not from the two-point function, but directly from the decay amplitudes (see 
also [34,35]). In addition, in Ref. [36] it was proposed to reconstruct the spectral density in the 
two-point function by using the maximum entropy method. This approach, in principle, also has 
the capability to address the problem of unstable states in lattice calculations. 

It turns out that the lowest-mass strongly interacting unstable particles in nature, the p(770), 
the A(1232)-resonance, etc 2 , due to their large width can not be identified with a clearly visible 
bifurcation in the energy levels. Namely, predicting the volume dependence of the pertinent 
energy levels by using the experimentally measured irn and vrA^ phase shifts, it is seen that the 
avoided level crossing is almost completely washed out. In this case, it is natural to ask, whether 
the resonance parameters which will be extracted by fitting Liischer's formula to the lattice data, 
can be determined at a reasonable accuracy and will be devoid of any bias. The issue of accuracy 
becomes particularly important since in forthcoming lattice calculations resonance parameters 
will be fitted to a few available data points. For example, in Ref. [37] the parameters of the 
p-meson have been determined performing a fit of Liischer's formula (in a moving frame) only to 
two data points. We believe that with more data on resonances expected to come, this question 
should be urgently addressed. 

In a previous paper [38] we have studied the problem in the case of the A-resonance. Invoking 
chiral effective field theory with explicit spin-3/2 degrees of freedom, we have parameterized the 
volume-dependent energy spectrum of the total Hamiltonian in terms of the A-resonance mass 
and width up to third order in the so-called small scale expansion (see, e.g. [39,40]). On the 
basis of a detailed analysis of the behavior of the two lowest energy levels, it was concluded that 
an accurate extraction of the A-resonance parameters is indeed a feasible task, despite the fact 
that the avoided level crossing is completely washed out. 

In the present paper we address the same problem within non-relativistic effective field theory 
(NR EFT) in a finite volume, which enables one to carry out the analysis in a more general 

2 We eschew here the ct(600) since at present there is no consensus about the precise nature of this resonance. 
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way 3 . The equation that determines the location of the eigenvalues of the Hamiltonian in this 
framework coincides with Liischer's formula. In order to facilitate the analysis, we further define 
a so-called probability distribution, which is constructed from the volume- dependent energies. 
The central observation is that the probability distribution in the vicinity of a resonance behaves 
much like the infinite-volume scattering cross section: it peaks at the resonance energy. The 
peak has approximately a Breit-Wigner shape, with the same width as the original resonance. 
We will show in the following that in case of a wide resonance, when the avoided level crossing 
is washed out, one still observes a clear resonance structure in the probability distribution after 
subtracting the background corresponding to the free motion of the decay products. This result 
unanimously supports the conclusion of Ref. [38]: the extraction of both the energy and width 
of the A-resonance from the volume-dependent spectrum by using Liischer's formula is feasible. 
Note also that, as shown in the present paper, this goal can be achieved even by fitting to the 
data for the lowest energy level alone. 

An important issue, which we do not discuss in the present paper, concerns the quark mass 
dependence of the resonance observables and the energy spectrum. For the values of the pion 
masses, which are used in present day calculations, the volume dependence of the energy levels 
may qualitatively differ from what happens at the physical value of the pion mass. For example, 
at higher pion masses the A-resonance becomes stable and the peak in the probability distribution 
degenerates into a 5-function. Our approach, combined with chiral perturbation theory, is flexible 
enough to describe this continuous transition. We, however, relegate a detailed discussion of this 
question to a separate publication [42]. 

The layout of the present paper is the following. In section 2 we consider the NR EFT for 
the 7iN system and derive Liischer's formula for spin-0 particle scattering on a spin-1/2 particle. 
Section 3 contains the reduction of Liischer's formula, using the cubic symmetry of the lattice. 
In section 4 we introduce the notion of the probability distribution and discuss its properties, 
as well as the infinite- volume limit. Finally, section 5 contains an analysis of (synthetic) lattice 
data, performed with the use of the probability distribution technique. We end with a summary 
and conclusions in sec. 6. Some technicalities are relegated to the appendices. 

2 Non-relativistic EFT for the pion-nucleon system in a 
finite volume 

Liischer's formula, which relates the infinite-volume elastic phase shift to the finite-volume two- 
particle energy spectrum, is derived in large volumes. Namely, the size of the three-dimensional 
box L, in which the two-particle system is placed, should be much larger than the typical scale 
M~ l set by the mass of the lightest particle (the pion in our case), in order to be able to discard 
all exponentially suppressed contributions at large L. For such large volumes, NR EFT provides 
an adequate description of the system at low energies. Liischer's formula within NR EFT can be 
straightforwardly obtained (see, e.g. [28]). In this section we briefly describe the generalization of 
the method to the case of particles with spin. Although the approach is completely general, below 
we shall focus on the example of pion-nucleon scattering. We shall use the covariant formulation 
of the NR EFT, introduced in Refs. [43]. The relativistic kinematics is taken into account 

3 Some of the results of this paper have been reported previously [41]. 
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automatically in this formulation that, in particular, may prove advantageous for generalizing 
Liischer's approach to moving frames. A recent general introduction to the NR EFT can be 
found, e.g. in Ref. [44]. 

The derivation consists of two parts. In the first part we set up the framework by considering 
the pion-nucleon scattering process in dimensionally regularized NR EFT in the infinite volume. 
To ease the notation, we suppress the isospin indices everywhere in the following, considering the 
scattering process in a channel with fixed total isospin. The non-relativistic Lagrangian takes 
the form 

C = $ f 2W 7T (id t - W n )$ + & 2W N (id t - W N )V + d , (1) 

where $ and ^/ denote the non-relativistic pion and nucleon fields, respectively. Further, W n = 
(M 2 — V 2 ) 1 / 2 and Wn = {m^ — V 2 ) 1 / 2 , with M w and m N the physical masses of the pion and the 
nucleon, in order. The ttN interaction Lagrangian Cj contains a tower of local 4-particle operators 
with increasing powers of space derivatives. The number of heavy particles is conserved. The 
coupling constants in £j encode the whole information about the high-energy behavior of the 
theory and are determined through matching to the effective-range expansion of the physical 
amplitudes. 

The non-relativistic pion and nucleon propagators are given by 

S M = o 7 \ — , \ 1 n ^ ' Sn ^ = o ~( \ ! \ 1 o ^ ' ( 2 ) 

2w w {p) ^(p) - p° - %0 2w N {p) w N {p) - p° - tO 

where w n (p) = (M 2 + p 2 ) 1 / 2 and wjv(p) = {m 2 N + p 2 ) 1 / 2 . 

The Lagrangian in Eq. (1) generates loops through the usual Feynman diagrammatic tech- 
nique. In order to ensure power counting and relativistic covariance, Feynman rules are sup- 
plemented with an additional prescription [43]: the integrands in all Feynman integrals are 
expanded in the inverse powers of masses, integrated by using dimensional regularization and 
finally summed up again to all orders. In the two-particle sector, which is considered here, this 
procedure is straightforward, since all loop contributions can be expressed through the basic 
bubble integral 



J{s) 



d u l 1 1 



(2tt) d 2w n (\)2w N (P - 1) K(l) - l°)(w N (P - 1) - P° + 1°) 



^( s ) mi j o\ A 1/2 (s,M 2 ,m 
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where A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx denotes the triangle function, D is the number 
of space-time dimensions and d = D — 1. 

Using canonical formalism, the full Hamiltonian of the nN system H = Ho + H/ can be con- 
structed from the Lagrangian (1). The scattering matrix T(z) is defined through the Lippmann- 
Schwinger (LS) equation 

T(z) = (-Hj) + (-H 7 )(-GoW)T(z) , (4) 

where Gq(z) = (z — Ho) -1 is the free resolvent of the irN system. Note that we have chosen to 
introduce negative signs in the above equation, in order to take into account the different sign 
conventions for the T-matrix in the field theory and in the potential scattering theory. 
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Next, we define the center-of-mass (CM) and relative momenta of the pion-nucleon pair, P 
and k, respectively, 

The pion-nucleon states are given by 

\(p N v),p w ) = \P,k,u), (6) 
where the index v labels the nucleon spin. The normalization of the states is fixed by 

(P', k', z/|P, k, v) = 2w N (p N )2w 7r (p 7r )5 v ^) d S d ('P' ' P)(2n) d 6 d (k' - k) . (7) 
We remove the CM momentum in the matrix elements by defining 

]d-r>> 



k;z) = J ^(P',k',v'\T(z)\0,k,v) 
h v , v {k!,k) = -J^(P',k!y\U\0,k,u). 



The LS equation in the CM frame takes the form 

Jd u" 



J {2n) d 2uj N {k') 2uj 7r (k") uj N {k") + uv(k") - z 

Since h u > u (k', k) is (an infinite) polynomial in momenta, by using Eq. (3) the above equation 
simplifies to 

t u , u (k', k; z) = h u , u (k', k) + 2 / dfi k " ^w(k', k") t u » u (k" , k; z) , (10) 

where <if2k» stands for the integral over the 3-dimensional solid angle and 

- k" \^\z\Mlm\) 

k _ |k"| 2z ■ { ] 

Note that, since the pion-nucleon loop in Eq. (3) is finite at d — > 3, one may set d = 3 in Eq. (10). 
All terms that vanish in dimensional regularization after performing the expansion in the inverse 
powers of masses are disregarded. The partial-wave expansion proceeds then in the standard 
manner. It is carried out in terms of spinor spherical harmonics defined by 

y L A(k,v) = ^(Lm^m s |JM>y^(k) Xm >) , (12) 



mim s 



where k = k/|k| and Yj m (k) and Xm a iy) are the Legendre spherical function and the two- 
component nucleon spinor, respectively. The quantity (Lmi \m s \JM) stands for the pertinent 
Clebsch-Gordan coefficient. 
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Introducing for convenience the projectors 



n^(k', k) = v'))*yjkk v) , (is) 



the partial wave expansion can be written as 

Wk',M = 4vr^n^(k',k)t JL (fc / ,A ; ; 



JL 



^v(k',k) = 47r^n^(k',k)^ L (fc / ,fc) . (14) 

JL 

On the mass shell, k = k' and z(k) = (M 2 + k 2 ) 1 / 2 + (m^ + &; 2 ) 1 / 2 , the scattering amplitude and 
the matrix element of the Hamiltonian can be expressed through the elastic scattering phases 
fijL,{k) i n a standard manner 

*„(*,*„(*)) = ^ eip|M f""V 

h JL (k,k) = 5^ tani Jt (*) . (15) 
k 

At the next step we consider the same system placed in a finite cubic box LxLxL. The Feynman 
rules in a finite volume remain the same, except that the momentum integration everywhere is 
now replaced by a discrete sum 

d d k" 1 . „ 27rn 



1 , » 27rn ~ , . 

^]T, k" = -,n6Z 3 . (16) 



Our aim is to find the finite-volume energy spectrum of the system described by the Lagran- 
gian (1). To this end, note that the location of the eigenvalues of the Hamiltonian in a finite 
volume coincides with the (real) poles of the operator T(z), defined by Eq. (4), in a complex 
,2-plane. After removing the CM momentum, this equation becomes similar to Eq. (9), but with 
the integration replaced through the momentum sum, as in Eq. (16). The ultraviolet divergence 
can be most conveniently tamed by analytic regularization [19]. In the following, we do not 
indicate the regularization explicitly. 

The resulting equation can again be expanded in partial waves. Since the rotational symmetry 
in the infinite volume is now broken down to a cubic symmetry, the partial wave expansion of the 
matrix elements of the operator T(z) will not be diagonal in J,L and M anymore. In order to 
ease the notations, it is useful to introduce the multi-index A = (J, L, M). Defining the operators 

nl A (k',k) = (yjX(u,"')yyj'kkv), (ir) 

the partial wave expansion can be written as 

Wk',k;z) = 4TrJ2^i A (k',k)t A , A (k',k;z) , 

A' A 



^v(k',k) = 4nJ2nt: u A (k',k)h A , A (k',k) . 



A'A 
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Further, since h A i A (k', k) is calculated from the non-relativistic Lagrangian at tree level, it coin- 
cides with its infinite- volume counterpart and is diagonal 

h A > A {k', k) = 5 A , A h A (k', k) = h JL (k', k)8j>j8 L , L 8 M , M . (19) 

Next, in analogy with the infinite- volume case, one may derive the equation for the on-shell 
quantities. To this end, note that the off-shell contribution to the LS equation is exponentially 
suppressed by the box size L, since the singular energy denominator is canceled in this contribu- 
tion (for the proof of this statement, see, e.g. [26]). Thus, up to these exponentially suppressed 
terms, the partial-wave expanded LS equation in a finite volume can be rewritten as 

k \ 

t A > A (k, k;z(k)) - 5 A > A h A (k,k) = - — -—} h A/ (k,k)M A > A >>{k)t A » A (k,k]z(k)), 

8nz(k) 

(20) 

where we have defined 

M A > A {k) - Mj>L>M>,JLM{k) - — |— Jj2_^2^ k „2 _ k 2 • ( 21 > 

k" v 

Note that in deriving Eqs. (20,21) we have brought the energy denominator of Eq. (9) to the 
non-relativistic form (q 2 — A; 2 ) -1 , eliminating the square roots by multiplying the numerator and 
the denominator in this equation by the same algebraic expression and neglecting off-shell terms, 
which are exponentially suppressed. 

The quantity A4j'L'M',JLM{k) is related to its counterpart for spin-zero particles [19], accord- 
ing to 

M J'L'M',JLM (k) = M L 'm',Lm(k) (L'm! -a\J'M')(Lm-a\JM) , (22) 

m'ma 

where [19] 

(— ) i ' L+L ' J p 

71 j=\L-L'\ s=-j K 



with k = kL/(2ir) and 



^ 2 ) = £t^# e*) 

nez 3 1 ; 



The coefficients CL> m ',js,Lm are expressed through the Clebsch-Gordan coefficients 



C Vml , js , Lm = lL '~ j+L \j i2L '^ ) ^ 3 1) + ^ (L'0j0\L0)(L'm'js\Lm) . (25) 
Due to Eq. (22) and the symmetry properties of M.L'm',Lm (see [19]), one finds that 

M-JLM,J'L'M' = M-J'L'M'.JLM ■ (26) 
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The quantity £a'a(&, k; z{k)) denned by the finite-volume LS equation (20) develops poles at 
the momenta where the determinant of this linear system of equations vanishes. Expressing 
hji(k,k) through the infinite-volume phase shift according to Eq. (15), we finally obtain the 
Liischer formula for pion-nucleon scattering 

det [tan 6 jii/(k)Mj>i/M',jLM{k) ~ ^j'j^ul^m'm] = . (27) 

This formula relates the location of the energy eigenvalues of the pion-nucleon system, placed in 
a finite box, to the infinite volume partial-wave phase shifts. 

In analogy to Ref. [19], it is possible to use the cubic symmetry on the lattice in order to 
achieve the partial block-diagonalization of the matrix A4j>L'M',jLM{k)- Such a reduction will be 
considered in the next section. 



3 Reduction of Liischer's formula 

In the infinite volume, the basis vectors of the irreducible representation D J of the rotation group, 
corresponding to the total momentum J, are given by \JLM) = ^2 ms \Lm^s)(Lrn^s\JM). 
Here, J = |,|,...,M = — J, ... J and L = J±|. The vectors \JLM) are also parity eigenvectors 
with the eigenvalue P = (— ) L . Below, we shall use a notation where the parity is explicitly 
indicated \JLM) = |JM) ± . 

In a finite volume, for the case of particles with the half-integer spin, the symmetry breaks 
down to 2 <g> 1S2, where 2 denotes the double cover of the cubic group containing 48 elements, 
no reflections included (see, e.g. [45]) and S 2 is the discrete group of space inversions. The 
irreducible representations of this group are Gf, Gf and H ± (see appendix A for details). The 
linear space spanned by the vectors | JM)^ 1 forms a basis of a reducible representation of the 
group 2 (8> 5*2. We denote the basis vectors, corresponding to the irreducible representations, as 

|r,a,J,n) ± , a = l,...dimr, n = 1, . . . N(T, J) . (28) 

Here, T = G±, G2 or H, N(T, J) denotes the multiplicity of the irreducible representation r in 
D J and the index a labels the vectors of a particular irreducible representation. 

The basis vectors \T,a, J^n)^ can be expressed through linear combinations of \ JM) ± 

\T,a,J,n) ± = J2^j2M\JM) ± . (29) 

M 

The matrix elements of the operator M.{k) in the new basis are given by 

± (r',a / , J',n'\M(k)\T,a, J,n) ± = ^^J'l'm'T c j2m Mj>v M >,jLM(k) . (30) 

M'M 

According to Schur's lemma, the operator A4(k) is partially diagonalized in the new basis 

± (T\a\J\n'\M(k)\T,a,J,n} ± = 8 rT 5 ala [Mi(k)}j n , Jlnl (31) 
and equation (27) is rewritten as 

ldet(taji6j fU (k)[Ml(k)} J'n'Jn - 5j>j5 n ' n = , (32) 
v=j>±\ r ^ / 
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r± 


J 


J' 


M T ± (k) 


Gf 


1/2 


1/2 


Woo 


Gf 


1/2 


7/2 


^f 1 w 40 


Gt 


7/2 


7/2 


Woo + ±fW 4 o + f>W 6 o 


Gf 


5/2 


5/2 


Woo-f W 40 


Gf 


5/2 


7/2 


± 60 7 f W 40 T 40 1 f w 60 


Gf 


7/2 


7/2 


Woo- ±f W 40 + f? W 60 


H± 


3/2 


3/2 


Woo 


H± 


3/2 


5/2 


T^W 40 


H± 


5/2 


5/2 


W 00 + fW 40 


H± 


3/2 


7/2 


7 ^^O 


H± 


5/2 


7/2 


T ^W 40 T^W 60 


H± 


7/2 


7/2 


Woo + # W 40 - | W 60 



Table 1: Non-vanishing matrix elements [.M±(/c)] j n ,j'n' for J, J' < 9/2 and n — n' — 1. The 
matrix is symmetric under J'n' <-» Jn. 

where the +/— sign in Ai± corresponds to even/odd V . 

In table 1 we list the matrix elements [-A/l^(/c)] j> n ',jn for J', J < |. Since the multiplicity 
N(T, J) = 1 for J < §, the indices n', n can be omitted in this table. The entries of the table are 
expressed through the following quantities [19] 

W m = ( 7 r 3 / 2 (2/ + l) 1 /V +1 )- 1 Z im (l;/, 2 ). (33) 
The construction of the basis vectors in case of arbitrary J is considered in appendix A. 

4 Probability distribution 

As mentioned in the introduction, in the vicinity of a narrow resonance the finite-volume energy 
levels of a two-particle system exhibit the peculiar behavior known as the avoided level crossing. 
Such a behavior, which is predicted by Liischer's formula, is schematically shown in Fig. 1. In 
this figure, we plot the relative momentum p, which is related to the CM energy E as E = 
{m 2 N +p 2 ) 1 ^ 2 + (M 2 +P 2 ) 1 / 2 , vs the box size L. The plateaus correspond to the resonance energy 
and the resonance width is determined by the minimal distance between the curves. 

It was, however, also mentioned above that for most physically interesting strong resonances 
the avoided level crossing is almost completely washed out from the spectrum due to the large 
width of a resonance. In this section, we describe a method that can be used to visualize the 
extraction of the resonance parameters from the two-particle spectrum even in this case. What 
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Figure 1: Schematic plot describing the construction of the probability distribution from the 
energy levels. 

is important is that the method does not contain any prior theoretical bias (e.g. does not use 
the resonance parameterization of the infinite- volume scattering phase as an input). 

Assume now that the volume-dependent two-particle spectrum is measured on the lattice. 
The probability distribution W(p) is constructed according to the following prescriptions: 

i) Choose the first N energy levels (e.g., N = 2 in Fig. 1); choose the interval L G [L\, L 2 ] and 
slice this interval into equal parts with length AL. For each value of L = L±, L = L + AL, 
etc determine p n (L), n = 1 ■ • • N. 

ii) Choose a corresponding momentum interval p e [pi,P2] an d introduce equal-size momentum 
bins with length Ap. 

iii) Count, how many times the eigenvalue p n {L), n = 1 • • • N is contained in a particular bin, if 
L runs from L\ to L 2 . This number gives the unnormalized probability distribution for the 
momentum bin chosen. Normalizing this distribution in the interval [pi,P2] yields finally 
the probability distribution W(p) we are looking for. The normalization condition is given 
by ' 

M P 

V W( Pl + kAp) Ap = 1 , M P = P2 ~ Pl - 1 . (34) 

Ap 

k=0 ^ 

It is clear that in case of a pronounced avoided level crossing, the probability distribution 
W(p) must be strongly peaked around the resonance energy. The exact shape of W(p) can be 
predicted on the basis of Luscher's formula. To this end, note that, in the limit of infinitesimally 
small L— and p— bins, W(p) is given by 

N 1 

n=l P ^ L > 
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where C is a normalization constant. 

For simplicity, we consider here the scattering in the partial wave with L = 1, J = | and 
neglect the (small) mixing to higher partial waves. Using table 1 and Eq. (32), the relation 
between the finite- volume energy spectrum and the phase shift S(p) takes the form 4 

S(p) = -0(k) + vrn, 0(«) = - arctan — — — — , k = — , (36) 

Z 00 (1;k 2 ) 2tt 

where the integer n labels the energy levels p n {L), which are the solutions of the above equation. 
Differentiating now Eq. (36) with respect to L and substituting into Eq. (35), we obtain 



71=1 ^ 



'L n (p) , 2tt5'(p) 



p P<j>'(K n (p)) 



(37) 



where K n (p) and L n (p) are the solutions of Eq. (36) for a given p. It is seen that W(p) defined by 
Eq. (37) is closely related to the so-called "density of states in a finite volume," see e.g. Ref. [46]. 

In the vicinity of the resonance, 8'(p) is strongly peaked. Substituting a Breit-Wigner param- 
eterization for S(p) and assuming that all other factors smoothly depend on the momentum p, 
we may verify that in the vicinity of the resonance the function W(p) follows the Breit-Wigner 
form for the scattering cross section, with the same width 5 . 

A useful parameterization of W(p) can be obtained by using the following approximation of 
the function 0(k), which is valid in a large interval of arguments [20] 

4>(k) = ttck 2 , c ~ 1 . (38) 

Solving Liischer's equation, we obtain 

L n (p) = 1 v/47r(7m - 5{p)) (39) 

and 

= £ f fy/S&EM + ■ MM ) . (40) 

In order to suppress the (large) background, related to the free motion of the 7rN pair, we consider 
in the following the so-called subtracted probability distribution W(p) — Wo(p), where Wo(p) is 
determined from Eq. (37) with 5(p) = and L n (p) corresponding to the free energy levels. 

It is interesting to consider the infinite- volume limit of the probability distribution. Note that 
in this limit the number of energy levels per fixed momentum bin goes to infinity. Consequently, 
the number of levels N in Eq. (36) can be chosen very large. In this case, in the expression for 
the quantity L n (p) which is determined through the solution of Liischer's equation 

L n (p) = —r 1 (™ - S(p)) , (41) 

p 



4 To ease notation, we do not attach indices L, J to this phase shift. 

5 Note that the location of the maximum in the probability distribution does not, in general, coincide either 
with the real part of the pole position in the amplitude, or with the solution of the equation 5{jpr) = ir/2. These 
three quantities agree only in the limit of the infinitely narrow resonance. 
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M L 

7t 



Figure 2: Volume dependence of the spectrum of the Hamiltonian (solid lines), predicted by 
using Luscher's formula with input experimental P33 phase shift. The center-of-mass momentum 
p vs the box size L is shown. For comparison, the free energy levels for the ttN system are given 
(dashed lines). The horizontal line marks the position of the A resonance. It is seen that the 
avoided level crossing is completely washed out. 

the phase shift obeys the inequality 5{p) <C im for the large majority of terms in the sum over 
energy levels. Expanding in a Taylor series, we get 

9tt 9tt 1 

L n ( P ) = -R n - -5{v)-j7-T + 0(5 2 ) , (42) 

with R n is the solution of the equation (f)(K n ) = nn. Using, in addition, n n = R n + 0(5), the 
unnormalized probability distribution for N — > 00 takes the form 

C-^^t^t^^f-m)^. (43, 

n=l n=l v ' x ' 

The first term exactly coincides with the free background. Subtracting this background and 
taking into account the fact that the quantity ^^^(^'(^n))" 1 does not depend on the phase 
6(jp), we obtain 

C-'Wip) - C^Woip) oc - (M - 5>(p)) . (44) 

p v v J 

In other words, in the infinite- volume limit this quantity is determined by the elastic phase shift 
alone. 



5 Analysis of synthetic data 

In this section we will implement the method discussed in the previous section for analyzing data. 
In the absence of lattice QCD data we will use synthetic data on the spectrum of the Hamiltonian, 
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Figure 3: Unsubtracted (left panel) and subtracted (right panel) probability distributions. Only 
the lowest energy level has been included in the analysis (N = 1). The solid lines correspond to 
the prediction made by using Liischer's formula with an approximation <P{k) ~ ttk 2 (see the text 
for details). A clear resonance-like structure is observed in the subtracted distribution. 

which are produced by using experimentally measured phase shifts [47] in Liischer's formula. If 
in the future unquenched lattice calculations are performed at the physical quark masses, the 
results must agree with the above synthetic data set. 

The calculated spectrum, obtained by substituting the resonant P 33 partial wave phase shift 
into Liischer's formula, is shown in Fig. 2. In this figure, the relative momentum of the 7riV 
system p(L), corresponding to the discrete energy levels of the Hamiltonian in a finite box, is 
plotted against the box size L (in units of M~ l ). It is seen that the structure of the energy levels 
is smooth: the avoided level crossing has been completely washed out due to the relatively large 
width of A. However, in the same figure we also plot the free energy levels, demonstrating that 
in the vicinity of the resonance a continuous rearrangement of the spectrum takes place. 

This rearrangement can be made explicitly visible by performing the analysis of the en- 
ergy spectrum, using the probability distribution method introduced above. Constructing, as 
described above, the unsubtracted probability distribution W(p) from the lowest energy level 
yields the plot shown on the left panel of Fig. 3. The resonance is seen as a barely distinguish- 
able shoulder around p ~ 0.22 GeV. This result, which obviously reflects the washing-out of the 
avoided level crossing in Fig. 2, casts justified doubts on the feasibility of a clean extraction of 
the A-resonance parameters from the data. The picture, however, completely changes once the 
subtraction of the background due to free irN pairs has been performed, see the right panel in 
Fig. 3. The resonance-like structure in the subtracted distribution is clearly visible, allowing one 
to finally conclude that the determination of the resonance parameters from the data is indeed 
possible. On both plots the solid curves correspond to the theoretical prediction made on the 
basis of Liischer's formula. In order to simplify the numerical calculations, the curves were con- 
structed by using the approximation <f)(n) ~ ttk, 2 , which works very well for all relevant values 
of the variable k. Since these curves are shown for the demonstrative purposes only, a better 
accuracy is not needed here. 

Up to now, we have dealt with the exact solution of Liischer's equation for the energy spectrum 
that on the lattice corresponds to measuring this spectrum at infinite accuracy and at all values 
of L. The situation in real calculations is different. Here one expects to get at most a few data 
points for different volumes. Our next aim is to mimic this situation in the calculations with 
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Figure 4: The same as in Fig. 3 but obtained only from 5 data points corresponding to M W L = 
1.4, 2.55, 3.7, 4.85, 6. The values of the function p(L) between the data points are obtained 
using Spline interpolation. Only the lowest energy level is analyzed (N = 1). The interval in 
the variable p is the same as in Fig. 3. It is seen that 5 data points do not provide a very good 
accuracy over this interval in L: the subtracted probability distribution is rather different from 
the theoretical prediction made on the basis of Liischer's formula. 

the synthetic data and check whether the extraction of the resonance parameters is still possible. 
The probability distribution method, which we are using, is equivalent to Liischer's approach 
and provides just a nice tool to visualize the final result. 

In order to achieve the goal formulated above, we first perform the analysis in the same 
momentum interval as in Fig. 3, but using 5 uniformly distributed data points, located at M n L = 
1.4, 2.55, 3.7, 4.85, 6. The values of p(L) between these data points were reconstructed by using 
the interpolation procedure with cubic splines. The result for the unsubtracted and subtracted 
probability distributions is shown in Fig. 4. It is clear that providing only 5 data points at this 
quite large interval of the variable L does not ensure a very high accuracy: the shape of the 
resonance comes out distorted. Moreover, one might expect that if even less data points are 
included in the analysis (see, e.g. [37]), the determined resonance parameters will include large 
systematic uncertainty which is very hard to control. 

There can be two possible ways out. One may try to gradually increase the number of data 
points, or one may try to reduce the size of the momentum and volume interval just to an 
immediate proximity of the resonance. Both possibilities have been tried, as shown in Fig. 5. 
The left panel of this figure corresponds to using 10 data points instead of 5 in the same interval, 
whereas the right panel corresponds to reducing the interval to M n L £ [1.9,4.5], respectively. 
In both cases one observes a clear improvement as compared to the case shown in Fig. 4. Note 
also that the data points should not not be uniformly distributed in L, but have indeed to be 
concentrated on both sides of the resonance. Otherwise, one may arrive to the picture shown 
in Fig. 6, where we show the probability distribution, obtained from the following data points: 
M n L = 1.4, 1.6, 1.8, 2.8, 3.9. It is seen that the probability distribution significantly deviates 
from the exact theoretical prediction in that part of the interval, where the data points are sparse. 

Finally, we have applied the method of probability distributions to the simultaneous analysis 
of the first two energy levels. The figure 7 contains the information about 10 data points, 
uniformly distributed in the interval M n L £ [2,6.5]. The resulting resonance shape is again in 
good agreement with the theoretical prediction, made on the basis of Liischer's formula. Note 
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Figure 5: The subtracted probability distribution for 10 data points equally distributed in the 
interval M n L = 1.4- • -6 (left panel). The same with 5 data points, M n L = 1.9 •• -4.5 (right 
panel). In both cases, only the lowest energy level has been included in the analysis (N = 1). 
The agreement with the theoretical curve, based on Luscher's formula, is much better than for 
Fig. 4 (right panel). 




p [GeV] 



Figure 6: A non-uniform choice of the data points M n L = 1.4, 1.6, 1.8, 2.8, 3.9. The momentum 
interval is chosen to be p £ [0.21, 0.35] GeV. It is immediately seen that the subtracted probability 
distribution deviates from the theoretically predicted behavior in the interval where the data 
points are sparse (large volumes or small momenta). 

that the resonance parameters extracted from the analysis of the different energy levels must of 
course coincide. In case the data from the excited levels are also available, checking the stability 
of the resonance parameters might enable one to verify a posteriori, whether the volumes, used 
in the calculation, are large enough to justify the application of Luscher's approach. 

Last but not least, the lattice data on the measured spectrum always come with errors. 
Our method provides an easy tool to render the analysis transparent in this case as well. The 
procedure is described below. 

We start from the data equidistantly distributed in the interval M n L £ [1.9,4.5] (5 or 10 
data points, as in Fig. 8). From each data point E{ = E(Lj) we further produce a statistical 
sample of 50 data points at a same Li, which are normally distributed around the central value E{ 
with the standard deviation o~i = eEj,. The figure 8 shows the probability distributions obtained 
from these randomly produced data, for three different values of e = 0.1%, 0.25%, 0.5%. Spline 
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p [GeV] 

Figure 7: The result of the analysis of the lowest two energy levels (N = 2) with 10 data points, 
M n L = 2 - ■ -6.5. The resonance structure reasonably reproduces the theoretical prediction on 
the basis of Liischer's formula. 

interpolation is used between data points. It is evident that the nice resonance structure, which 
was seen in the probability distribution (data with no errors), is washed out already at quite 
small values of the relative error assigned. Interestingly enough, the increase of the number of 
data points does not lead to an improved accuracy. The reason for this is that we treat the 
neighboring data points to be statistically independent. If the distance between the neighboring 
points is decreased, the fluctuations in the derivative of the eigenvalues increase and this leads 
to the increase of the statistical noise in the probability distribution (see Fig. 8). 

These statistical fluctuations can be suppressed to some extent if we perform a smooth (e.g. 
polynomial) interpolation of random data prior to calculation of the probability distribution. 
This is demonstrated in Fig. 9 where the random data on the spectrum were first fitted as 
p(L) = a Q + a\L + (I2L 2 + a 3 L 3 prior to producing the probability distribution. The situation 
clearly improves, especially in the case of 10 data points. Still, from Figs. 8 and 9 one has to 
conclude that a very accurate measurement of the spectrum is indeed needed to reliably extract 
the properties of the A-resonance. 

6 Conclusions 

i) Within the covariant non-relativistic effective field theory [43] we have derived Liischer's 
formula for scattering of spin- 1/2 and spin-0 particles. The partial- wave expansion is per- 
formed, and the cubic symmetry on the lattice is used to reduce the resulting matrix equa- 
tions. 

ii) The notion of probability distribution for the finite-volume spectrum of the Hamiltonian is 
introduced. It is shown that near the resonance energy the probability distribution behaves 
similar to the scattering cross section in the infinite volume: it produces a Breit-Wigner 
peak at the resonance energy with the same width. 

iii) The probability distribution, which is directly constructed from the energy levels, does not 
contain any prior bias. For this reason, the analysis carried out with the use of this method 
can be used to judge whether a clean extraction of the resonance parameters from the 
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Figure 8: Probability distributions, obtained from the data that contain errors (see the text for 
the details). The central data points are the same as in Fig. 5. As seen from the figure, the 
resonance structure is effectively washed out already at a relative error of 0.5% in the data or 
even earlier. 



available data is possible. 

iv) The probability distribution does not carry more or less physical information than Luscher's 
formula. The advantage of the probability distribution is its visual transparency. The choice 
of the method in the actual analysis is dictated by convenience. 

v) In the present paper we apply the method of probability distribution to the case of the A- 
resonance. We observe that the distribution after subtracting the background corresponding 
to the free motion of the ttN pair develops a nice resonance structure in accordance with 
the exact prediction based on Luscher's formula - even though the avoided level crossing is 
completely washed out. 

vi) It is possible to achieve a satisfactory description of the resonance position and shape even 
with few data points, provided they are chosen close enough to the resonance and are 
measured very accurately. Measurement of only the ground state suffices (inclusion of the 
excited levels provides additional check on the results). To conclude, the results of the paper 
clearly demonstrate that the extraction of the resonance parameters from the measurement 
of the finite- volume energy spectrum by using Luscher's method is indeed a feasible although 
difficult task. 
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Figure 9: The same as in Fig. 8 at a relative error e = 0.5%. A polynomial fit to the random 
data was used prior to producing the probability distribution. The improvement as compared to 
Fig. 8 is clearly visible. The peak is washed out approximately at e ~ 0.75% ■■■!%. 
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A Basis vectors of the cubic group 

In this appendix we collect the basic formulae for the representations of the cubic group, which 
are needed for partial diagonalization of Luscher's formula. 

Consider first shortly the case without spin. The group O of the symmetries of the 3- 
dimensional cube has 24 elements (no reflections included), which fall into 5 conjugacy classes: 
/ (identity), 36*2, 8C3, 6C4 and 6C' 2 (see, e.g. [45,48]). We find it convenient to parametrize 
the group elements by specifying three Euler angles a,/3, 7. Alternatively, the element can be 
parametrized, e.g., by specifying the axis n and the rotation angle u. The table A.l collects the 
values of the group parameters. 

The irreducible representations of the rotation group D L , L = 0, l,--- are defined in the 
2L + 1-dimensional space spanned on the basis vectors \LM). These representations are reducible 
under the cubic group O and can be decomposed into the irreducible representations of the latter 
denoted by T = Ai, A 2 , E, T\ and T 2 . The dimension of these representations N(T) is equal to 
1,1,2,3,3, respectively. 

In order to construct the basis of the irreducible representations, we consider the linear 
operator P T a f ■, whose matrix elements in the space spanned by the vectors \LM) are given by 

24 

{Plt)MM> = J2( R iU D MM>(ai,Pi,li) , (A.l) 
i=l 

where Dj^ M ,(ai, A;7i) are Wigner D-functions, and (Rj) a /3, a,P — 1- ■ ■ N(T) denote the matri- 
ces of the irreducible representations of cubic group 

A\\ A\ is the trivial 1-dimensional representation E4 — 1. 
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A2. Ri = — 1 for the conjugacy classes 6C4 and 6C 2 , Ri = 1 otherwise. 
E: The matrices in this representation are two-dimensional and real: 



Ri — 


1 for 


i = 1,22,23, 


24, 




JXi — 


cr 3 for 




8 1 Q 




Ri = 


71 

— cos — 
3 


71 

1 + z sin - cr 2 




for 


f = 2,5,6,9, 


Ri = 


7T 

-cos- 


7T 

1 — i sin — CT9 

3 


for 


2 = 3,4, 7,8, 


Ri = 


71 

-cos- 


71 

0-3 - sin - 0-1 


for 


i = 10,11,16,17, 


Ri = 


7T 

— cos — 
3 


7T 

0-3 + sin - cti 


for 


i = 12, 13,20,21. 



(A.2) 

T\\ (Ri) a /3 = exp( —in®Ju)i ) = cosoj i 8 ,^-\-{l — cos Wj)tia - sin ujie a ^n^\ where (J 7 ) a /3 = 

V /a/3 

—ie a f3-y denote the group generators. 

T 2 : The matrices are the same as in the irreducible representation T 1; except the change of sign 
for the conjugacy classes 6C4 and 6C 2 . 

The basis vectors of the irreducible representations are obtained by acting with the linear 
operator given by Eq. (A.l) at a fixed (3 and varying a on an arbitrary vector <pM from the space 
spanned by the vectors \LM) 

L 

{e r a L ' p ) M =M Yl ( P ap L )MM'<j>M> , a = 1 • • • N T , (3 fixed , (A.3) 

M>=~L 

where M denotes the normalization constant. It is fixed so that the basis vectors obey the 
orthonormality condition 

Ys^Tm (£ L ' P )m = Wrr ■ (A.4) 

M 

If the representation T is not contained in D L , the action of the projection operator on 0m gives 
0. The equations (A.3) and (A.4) do not fix a common phase of the basis vectors, belonging to 
the same representation labeled with L,T - this can be freely chosen. If a representation T enters 
more than once in D L , an additional orthogonalization of the basis vectors, belonging to the same 
representation, is necessary. Inclusion of parity is trivial. The basis vectors are simultaneously 
the eigenvectors of S2 with the eigenvalue P = (— ) L . 

In table A.2 we list the basis vectors of the irreducible representations of the cubic group up to 
L — 4, obtained from Eq. (A.3). The phases are chosen so that after the partial diagonalization 
of Luscher's equation, the entries of table E.2 in Ref. [19] are reproduced. 

Note also that our basis differs from the one given in Ref. [49] and can not be reduced to it 
with a single unitary transformation for all L. We have also checked that the use of the basis 
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from Ref. [49] does not lead to Eq. (31): e.g., the different irreducible representations turn out 
not to be orthogonal. 

Our basis vectors can be transformed into those listed in Ref. [50] by unitary transformations. 
Since in that article the basis vectors of the same irreducible representation, belonging to different 
values of L, are not fully displayed, we can not carry out this comparison to the end. 

In order to include the particles with spin, one has to consider the double cover of O denoted 
by 2 0, which can be constructed by adding a negative identity for ±2ir rotations to the group 
O [45]. One ends up with a group of 48 elements divided into 8 conjugacy classes (see table A. 3) 
and, accordingly, with 8 irreducible representations. In addition to the previously considered 
5, one has 3 new representations denotes as T = G\,G2,H with the dimension N(T) = 2,2,4, 
respectively. In case of the half-integer total momentum J, one needs to consider only these 
additional even-dimensional representations. 

The matrices of the irreducible representations Gi,G2,H are given by [45] 

G\\ {Ri) a p — expl — | nWffWi J = 5 a p cos ^ — in^cr^sm ^. 

G 2 : the matrices are the same except change the sign in the conjugacy classes 6Cg, 6C' 8 and 
12C 4 . 

H: the matrices {Ri) a p = expf — zn^ji iOi J where J** denote the group generators in spin- 

V / a/3 

3/2 case. 

The counterpart of Eq. (A.l) in case of the half-integer total momentum J is 

48 

{P V J)mm> = ^(R^U D J MMI {a l: Pi, 7i ) . (A.5) 
i=i 

Acting with the linear operator P a ^ on an arbitrary linear combination of the basis vectors 
| JM)^ defined in section 3, we obtain the basis of the irreducible representations G±, G 2 and H. 
Inclusion of parity is again trivial, since the basis vectors \ JM) ± are the eigenvectors of parity, 
with the eigenvalue P = ±1 = (— ) L . Up to the value J = |, these basis vectors are listed in 
table A.4 

Using this basis to partially diagonalize Liischer's equation, we finally arrive at the results 
displayed in table 1. 
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— 7T 


-tt/2 


tt/2 


-tt/2 




18 


(1,1,0) 


— 7T 


-tt/2 


— 7T 







19 


(1,-1,0) 


— 71 





TT 


-tt/2 




20 


(1,0,1) 


— 71 





tt/2 


— 7T 




21 


(-1,0,1) 


— TT 





-tt/2 


— 7T 


3C 2 


22 


(1,0,0) 


— TT 


7T 


7T 







23 


(0,1,0) 


— TT 





— 7T 







24 


(0,0,1) 


— 7T 








— 71 



Table A.l: Parameterization of the elements of cubic group. The vector n should be normalized 
to unity. 
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r 


L 


a 


Basis vectors 


At 





1 


|0,0) 


rr 


1 


1 


^(|1,-1>-|M» 






2 


^(|i,-i) + |M)) 






3 


|i, o> 




2 


1 


_-^(|2,-l) + |2,l» 






2 


^=(|2,-1)-|2,1» 






3 


-i(|2,-2)-|2,2» 


E + 


2 


1 


|2,0) 






2 


^(|2,-2) + |2,2» 


Tr 


3 


1 


f (|3,-3)-|3,3»-^ (|3,-1>-|3,1» 






2 


=^(|3 J -3> + |3 ) 3»-^(|3,-l> + |3,l» 






3 


|3,0) 




3 


1 


(|3,-3)-|3,3})-f (|3,-1>-|3,1» 






2 


•AT ■ /F 

=^(|3 ) -3> + |3,3» + ^(|3 I -l) + |3,l» 






3 


^(|3,-2> + |3,2» 


A~ 2 


3 


1 


^(|3,-2)-|3,2)) 


7f 


4 


1 


-i(|4,-3) + |4,3))-f (|4,-1) + |4,1)) 






2 


i(|4,-3)-|4,3»-^(|4,-l)-|4,l)) 






3 


^(|4,-4)-|4,4)) 




4 


1 


^(|4,-3> + |4,3»-±(|4,-l> + |4,l» 






2 


^(|4,-3}-|4 ) 3)) + i(|4,-l)-|4,l)) 






3 


^(|4,-2)-|4,2)) 


E + 


4 


1 


-jf (|4,-4) + |4,4» + ^|4,0) 






2 


-i(|4,-2) + |4,2)) 


At 


4 


1 


(|4,-4> + |4,4)) + *fi|4,0) 



Table A. 2: The normalized basis of the irreducible representations of the cubic group: integer 
values of the angular momentum. 
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Class 


% 


n 


UJ 


Class 


i 


n 


UJ 


/ 


1 


any 





8C 3 


28 


(1,1,1) 


4tt/3 


6C 4 


2 


(1,0,0) 


TT 




29 


(-1,1,1) 


4tt/3 




3 


(0,1,0) 


TT 




30 


(-1,-1,1) 


4tt/3 




4 


(0,0,1) 


TT 




31 


(1,-1,1) 


4tt/3 




5 


(1,0,0) 


— TT 




32 


(1,1,1) 


-4tt/3 




6 


(0,1,0) 


— TT 




33 


(-1,1,1) 


-4tt/3 




7 


(0,0,1) 


— TT 




34 


(-1,-1,1) 


-4vr/3 


6C 8 


8 


(1,0,0) 


tt/2 




35 


(1,-1,1) 


-4tt/3 




9 


(0,1,0) 


tt/2 


12C 4 


36 


(0,1,1) 


TT 




10 


(0,0,1) 


tt/2 




37 


(0,-1,1) 


TT 




11 


(1,0,0) 


-tt/2 




38 


(1,1,0) 


TT 




12 


(0,1,0) 


-tt/2 




39 


(1,-1,0) 


TT 




13 


(0,0,1) 


-tt/2 




10 


(1,0,1) 


TT 


6C 8 


14 


(1,0,0) 


3tt/2 




41 


(-1,0,1) 


TT 




15 


(0,1,0) 


3tt/2 




12 


(0,1,1) 


— TT 




16 


(0,0,1) 


3tt/2 




13 


(0,-1,1) 


— TT 




17 


(1,0,0) 


-3tt/2 




44 


(1,1,0) 


— TT 




18 


(0,1,0) 


-3tt/2 




45 


(1,-1,0) 


— TT 




19 


(0,0,1) 


-3tt/2 




16 


(1,0,1) 


— TT 


8C 6 


20 


(1, 1, 1) 


2tt/3 




47 


(-1,0,1) 


— TT 




21 


(-1,1,1) 


2tt/3 


J 


48 


any 


2tt 




22 


(-1,-1,1) 


2tt/3 












23 


(1,-1,1) 


2tt/3 












24 


(1,1,1) 


-2vr/3 












25 


(-1,1,1) 


-2tt/3 












26 


(-1,-1,1) 


-2tt/3 












27 


(1,-1,1) 


-2tt/3 











Table A. 3: Parameterization of the elements of the double cover of the cubic group. The vector 
n should be normalized to unity. 
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j 


a 


Basis vectors 




1/2 


1 

2 




II I\± 
1 2 2/ 

II _I\± 
1 2 21 


G f l ± 


7/2 


1 


VIE i7 

6 1 2 


_ 7\± , yJT |7 1\± 
2' 6 1 2 2' 






2 


-%/21 1 7 


1\± VV> 1 7 7\± 
2' 6 1 2 2' 








6 1 2 


G2± 


5/2 


1 


Y.30 1 5 
6 1 2 


_ 3\± \/6 | 5 5\± 

2 / 6 1 2 2 ' 






2 


V6 1 5 
6 1 2 


_ 5\± _L vJO 1 5 3\± 
2 ' 6 12 2' 


G2 ± 


7/2 


1 


1 1 7 

2 1 2 


3\ ± _i_ \/3 1 7 5\ ± 
2' 2 1 2 2' 






2 


_yj |I 

2 1 2 


_ 5\± _|_ 1 I 7 3\± 
2 / 2 1 2 2 ' 




3/2 


1 

2 
3 
4 




1- 

1 2 21 
1- 

1 2 2 ' 

1- — -) ± 

1 2 2' 

13 3\± 
1 2 2' 


#± 


5/2 


1 

2 


-v^ 1 5 
6 1 2 


5\± \/6 I 5 3\± 
21 6 1 2 2' 

15 1\± 
1 2 2' 






3 




1 5 1\± 
1 2 2' 






4 


a/6 I 5 


3\± , \/30 15 5\± 






6 1 2 


21 ^ 6 1 2 2' 


if* 


7/2 


1 


1 1 7 

2 1 2 


5\± i \/3 1 7 3\± 
2' ~ r 2 1 2 2' 






2 


v^T 1 7 


7\± \/T5 1 7 1\± 
2/ 6 1 2 2/ 








6 1 2 






3 


\/l5 1 7 


1\± 1 \/2T 1 7 7\± 

2/ 6 1 2 2/ 








6 1 2 






4 


V3 1 7 


3\± , 1 1 7 5\± 






2 1 2 


21 ' 2 1 2 21 



Table A. 4: Basis functions of irreducible representations of 2 in terms of the functions \JM) 
defined in section 3. 
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